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We present a theoretical study of the dispersion relation of surface plasmon resonances of meso- 
scopic metal-dielectric-metal microspheres. By analyzing the solutions to Maxwell's equations, we 
obtain a simple geometric condition for which the system exhibits a band of surface plasmon modes 
whose resonant frequencies are weakly dependent on the multipole number. Using a modified Mie 
calculation, we find that a large number of modes belonging to this flat- dispersion band can be ex- 
cited simultaneously by a plane wave, thus enhancing the absorption cross- section. We demonstrate 
that the enhanced absorption peak of the sphere is geometrically tunable over the entire visible 
range. 



In recent years there has been a growing interest in sur- 
face plasmon resonances (SPRs) of metal-dielectric struc- 
tures as means to concentrate electromagnetic (EM) field 
in subwavelength volumes. The geometric tunability 
and the enhanced field intensity at the metal-dielectric 
interface have led to a number of theoretical and ex- 
perimental demonstrations of SPR-assisted EM energy 
focusing^, Q, 0. In particular, extensive studies have 
been conducted on SPRs of nanoshells — thin metal shells 
surrounding sub-micron dielectric cores. Both single^ 
and concentric ^5] nanoshell systems have been addressed. 
It was shown that the resonant frequencies, where inci- 
dent EM radiation is efficiently focused to the nanoshells, 
are tunable from the visible to near-infrared. Moreover, 
the large EM field concentration at the metallic shell pro- 
duces a giant surface enhanced Raman scattering y. 

The size scale of a typical nanoshell is of the order of 
lOOnm or less, a small fraction of the excitation wave- 
length. As a result, only a single, low order multipole 
resonance plays a vital role in channelling the incident 
EM energy towards the metallic shell. On the other hand, 
mesoscopic layered particles open up the possibility for si- 
multaneous excitation of a large number of SPRs, further 
increasing the efficiency of the EM field focusing. This, 
in turn, enhances optical phenomena such as absorption 
and nonlinear response. 

In this Letter, we show how to utilize dispersion en- 
gineering to enhance the EM field focusing and the ab- 
sorption cross-section of mesoscopic multilayered spheres. 
The system consists of a spherical resonator comprised 
of thin, alternating layers of dielectric and metal shells 
around a concentric metal core. The composite particle 
is embedded in a homogeneous, isotropic dielectric host 
with permittivity eo- More specifically, we address metal- 
dielectric-metal (MDM) microspheres: metallic cores 
comparable in size to optical wavelengths, surrounded 
by one sequence of lossless dielectric shell of thickness L 
followed by a metal shell of thickness T. We show that it 
is possible to obtain a band of SPRs with nearly identical 
resonant frequencies (flat-dispersion band) by adjusting 
the above geometric parameters. By solving Maxwell's 



equations using a spherical multipole expansion, we cal- 
culate the dispersion relations for the fiat band and the 
resultant enhanced absorption cross-section. 

To analyze the SPRs excited in the system, we initially 
address the dispersion relations of a simplified geometry 
in which the outer metal shell is infinitely thick. This 
model is complementary to a metal nanoshell on a di- 
electric core 4J, the DMD. Unlike the DMD, the presence 
of a properly designed dielectric shell in the MDM is nec- 
essary and sufficient for achieving a flat-dispersion band. 

We first consider a metal characterized by a lossless 
Drude model, e^(a;) = 65 — cj^/cc;^, where 65 is the con- 
tribution of inter-band transitions and cOp is the plasma 
frequency. The dielectric shell has a real permittivity Cd- 
The eigenfrequencies of the simplified MDM system are 
determined by the eigenmode equation 
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where 77 = 1 for transverse electric (TE) modes and 
T] = Cd/cm for transverse magnetic (TM) modes, and 
km,d = s/^^i^ ^/c. The radius of the inner metal core 
is R and S = R -\- L. The spherical Bessel and Hankel 
functions of the first kind of integer order / are denoted 
respectively by ji and hi , and the prime denotes differen- 
tiation with respect to the argument. We limit our dis- 
cussion to large enough particles, satisfying 1 <C kdR and 
1 <C \km\R- Using a real e^(a;) leads to straightforward 
analytic solutions of Eq. (^) with real 00. We show later 
that using a realistic (lossy) metal only modifies quan- 
titative aspects of the flat band, while its fundamental 
physical origins remain unaltered. 

For each polarization there is an infinite number of 
solutions to Eq. (QJ, each characterized by a multipole 
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number, / and n > roots, the latter yielding a radial 
excitation number (i.e. band index). We first examine 
the resonant modes for asymptotic limits of Eq. For 
high-order multipoles satisfying / ^ l^ml'S' and I ^ k^S 
we expand ji and hi to obtain the resonance condition 



(2) 



for TM polarization. No similar resonance condition ex- 
ists for TE polarization. We next look for an expression 
for L such that Eq. ((2|) is also satisfied for the TM mode 
with the smallest multipole number. For / = 1 where 
expansions of ji and hi may be applied, Eq. (QJ reduces 
to 



(2n - l)7r = 2kdL, 
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with n > 1. 

The high order multipole modes of Eq. (|2j belong to 
the n = 1 radial excitation band, hence by setting n = 1 
in Eq. 0, we obtain the geometric condition 
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where Xgp = 27r c/cOgp and cOgp is given by emi^sp) -\- = 
0. The position of the flat band given by Eq. ^ and the 
geometric condition, Eq. do not depend on R and 
are identical to those for the planar MDM geometry 7J. 

Figure ^a) shows the resonant frequencies of the sys- 
tem plotted as function of the multipole number To 
approximate the optical response of silver, we have cho- 
sen 65 = 5.1 and hcOp = 9.1eV|8|]. The dielectric shell 
has Cd = 3.53, corresponding to the measured value of 
amorphous titania. For clarity, we have chosen a high 
dielectric constant to place uosp well away from the bulk 
plasma resonance at cOp. Using Eq. L"^ = 53.4nm, and 
we set R = 500nm. From the figure, it is clear that there 
is a band of TM modes whose frequencies are weakly de- 
pendent on the multipole number. This fiat-dispersion 
band is near ijJgpl^p = (^6 + ^d)~^^^ ^ 0.34 as expected, 
and the width of the band is given by 5ix)/uOsp = 0.014 
where Sco is the difference of the largest and the small- 
est frequencies in the band. For comparison. Fig. ^b) 
shows the dispersion relation of a dielectric sphere of ra- 
dius R-\- L embedded in a metallic host. 

By analyzing the dispersion relation of MDM spheres 
with various core radii, we have found that when R > 
lOOnm, Sco depends weakly on the core radius R. For ex- 
ample, Suo/uOsp = 0.018 for R = lOOnm and Suo/uOsp = 
0.014 for R = lOOOnm. We have observed that the 
material dispersion of the metallic medium affects the 
fiatness 6uj/uJsp more significantly; a larger and positive 
dem/duj\^^^ leads to a fiatter dispersion relation. 

Next, we analyze the absorption cross- section of the 
MDM sphere with a metal shell of finite thickness, and 
study the coupling between the incident plane wave and 
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FIG. 1: (a) Dispersion relations of TM (squares) and TE 
(crosses) modes of the MDM sphere. Near the plasmon fre- 
quency oosp (solid line), there is a flat band of TM modes 
weakly dependent on multipole number, (b) Dispersion rela- 
tion of same structure without metal core. The insets show 
schematic cross-sections of the systems. 



the flat-dispersion modes. Such spheres may be real- 
ized by chemical synthesis of layered metallodielectric 
particles [3 |^. We have developed an exact and nu- 
merically stable algorithm for calculating the absorption 
cross-section of multilayered spheres based on previous 
publications 0, O, ^3 Q- We also account for 
absorption losses in the metal by modifying the Drude 
model to Cmi^) = — + iVuo)~^ ^ with F describ- 

ing the electron relaxation rate. We set hV = 0.021eV 
and eo = 1. 

A sharp absorption peak near Xsp is seen in Fig. Elja). 
For a core oi R = 500nm this is maximized when L = 
62nm and T = 77nm, and it is positioned at A = 428nm. 
The deviations of the peak from Xgp = 401nm and of 
the optimal value of L from = 53.4 are explained as 
follows: In addition to SPRs of the flat-dispersion band, 
SPRs of the outermost metal-host interface contribute 
to this peak. In fact, the two SPR branches are coupled 
through the finite metal shell. These plasmon hybridiza- 
tions lead to slight modifications in peak position as well 
as in the optimal dielectric shell thickness. 

The degree of coupling between the incident plane wave 
and the fiat dispersion band is tuned by adjusting T, the 
thickness of the outer metal shell. As we show later, the 
value of L may be varied in the vicinity of to spec- 
trally tune the absorption maximum while maintaining 
a nearly fiat dispersion band. For each such value of L 
there exists a T which optimizes the coupling. Neverthe- 
less, if the metal shell is made too thin, strong plasmon 
hybridization will eventually distort the fiat band. 

The mode decomposition of the spectrum, shown in 
Fig. Eljb), confirms that a large number of multipoles 
share nearly identical resonant frequencies and shows 
that the first 13 multipoles are excited concurrently by 
the incident plane wave. Higher order multipoles have 
negligible contribution to the absorption peak. Com- 
pared to the absorption spectrum of a metal sphere and of 
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FIG. 2: (a) Optimized absorption cross-section spectrum of 
MDM sphere with R = 500nm. (b) Mode decomposition of 
the spectrum from / = 1 (bottom) to / = 15 (top), showing 
contributions from TE modes (blue hatches), and TM modes 
(red hatches), (c) Absorption spectrum of a metal sphere of 
radius R-\- L -\-T. (d) Absorption spectrum of a sphere with 
dielectric core of radius R-\-L and metal shell of thickness T. 
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FIG. 3: Absorption spectra of silver-titania MDM sphere 
with R = 500nm and various dielectric shell thickness L. The 
value of T is chosen to maximize the peak height. 

a core-shell of the same size, plotted in Fig. |2Jc) and (d), 
the absorption cross-section of the MDM near A = 428nm 
is enhanced by a factor of 10 and 4, respectively. The en- 
hancement factor relative to the maximal absorption is 
seen to be ~ 4 for both cases. The enhanced absorption 
exists for all core radii R ^ lOOnm. The value of R af- 
fects the peak height and the number of simultaneously 
excited SPRs. Other spectral characteristics, such as the 
peak position, its FWHM, and the optimal values of L 
and T, are insensitive to the core size. For example, a 
change of a few hundred nanometers in the core radius al- 
ters the latter by less than several nanometers. For radii 
smaller than lOOnm, the absorption broadens, splitting 
into multiple peaks. Below R = 50nm, the effect of the 
metal core on the absorption diminishes, the spectrum 



rapidly approaching that of a DMD nanoshell. 

The enhanced absorption is also demonstrated numer- 
ically using the experimentally obtained dielectric func- 
tion for silver [l^. For R = 500nm, the absorption cross- 
section peak is maximized for L = 94nm and T = 31nm, 
and it is at A = 533nm with 64nm FWHM as shown in 
Fig. |31 The absorption enhancement factors relative to a 
solid metallic sphere and to a core-shell sphere of a com- 
parable size are 7 and 3, respectively. While the spectral 
position of the flat dispersion (i.e. Usp) depends only on 
the optical properties of the constituent media, the ab- 
sorption peak position can be tuned by adjusting L. As 
discussed above, an arbitrary L does not produce a flat 
dispersion. However, the slope of the band, uoi — ui-i, is 
generally small as long as L ~ L"^. Thus, it is possible to 
vary the value of L by tens of nanometers and shift the 
resonant frequencies, and still retain a significant overlap 
of the absorption cross-section peaks of various multi- 
poles. Figure 3 shows the geometric tunability of the 
enhanced absorption over the visible range. 

In summary, we have analyzed the dispersion relation 
of mesoscale MDM spheres, and derived the condition for 
which the system exhibits a band of TM modes whose 
eigenfrequencies are weakly dependent on the multipole 
number. The numerically obtained absorption spectra 
exhibit enhanced absorption as a consequence of the si- 
multaneous excitation of a large number of SPRs. Since 
the enhancement implies a large EM field concentra- 
tion, these results suggest that the flat-dispersion MDM 
spheres may be used to enhance nonlinear optical phe- 
nomena, and more generally, they may lead to a new 
generation of mesoscale plasmonic systems in which nu- 
merous modes are excited to manipulate SPR-assisted 
light-matter interactions. 
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